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Abstract 

We present a complete systematics (excitation function, impact parameter, system size, isospin 
asymmetry and equations of state dependences) of global stopping and fragments production for 
heavy- ion reactions in the energy range between 50 and 1000 MeV/nucleon in the presence of 
symmetry energy and isospin-dependent cross-section. It is observed that the degree of stopping 
depends weakly on the symmetry energy and strongly on the isospin-dependent cross-section. 
On the other hand, the symmetry energy and isospin-dependent cross-section has an effect of 
the order of more than 10% on the emission of light charged particles (LCP's). It means that 
nuclear stopping and LCP's can be used as a tool to get the information of isospin-dependent 
cross-section. Interestingly, the LCP's emission in the presence of symmetry energy is found to be 
highly correlated with the global stopping. 
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I. INTRODUCTION 



One of the goals of heavy-ion collisions (HIC) at intermediate energies is to extend 
the knowledge of hot and dense nuclear matter to the extreme conditions. In the past, 
these studies were focused on multifragmentation, that constitutes fragments of all sizesjl]. 
Additional promising observable for the understanding of nuclear equation of state is the 
anisotropy in the momentum distribution that includes the directed in-plane flow (bounce 
off ) as well as out of plane flow (squeeze out) 2], ^. The absolute values of the flow 
results from the interplay between the attractive mean field and repulsive nucleon-nucleon 
scatterings. This interplay is also responsible for the transition from a fused state to 
total disassemble one. The another phenomena linked with the above interplay is the 
global stopping of nuclear matter. Recently, Puri and co-workers |4|, tried to correlate the 
multifragmentation with global nuclear stopping. Their findings revealed that light charged 
particles acts in a similar fashion like anisotropy ratio. They, however, did not take isospin 
of the system into account. 

Following the recent development of radioactive beam facilities in many parts of the world, 
it became possible to study the neutron (or proton) rich nuclear collisions at interme- 
diate energies. Therefore, for a meaningful investigation, one should include the isospin 
dependence of the field. As pointed out by Bauer jj], nuclear stopping at intermediate 
energies is determined by the mean field as well as by the in-medium NN cross-sections. 
Unfortunately, his calculations were silent about the symmetry potential. The recent work 
of many authors suggest that the degree of approaching isospin equilibration helps 

to probe the nuclear stopping in heavy-ion collisions. In Ref.!?], isospin dependence of 
cross-section was investigated in nuclear stopping. In a recent communication [s], authors 
studied the behavior of excitation function Q^z/nucleon and concluded that Qzz/nucleon 
can provide information about the isospin dependence in term of cross- sect ions. Several 
more studies have also focused in the recent years on the isospin degree of freedom 9|. 
We wish to focus on the systematic study of isospin dependence and will focus on the 
relation between light charged particles and equilibration of the reaction using isospin- 
dependent quantum molecular dynamics. 

This study is done within the framework of isospin-dependent quantum molecular dynamics 
model that is explained in section-II. The results are presented in section-Ill. We present 
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summary in section-IV. 



II. ISOSPIN-DEPENDENT QUANTUM MOLECULAR DYNAMICS (IQMD) 
MODEL 



The isospin-dependent quantum molecular dynamics (IQMD)[10| model treats different 



charge states of nucleons, deltas and pions explicitly 



111 ], as inherited from the VUU model 



12|. The IQMD model has been used success: 



observables from low to relativistic energies 
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uUy for the analysis of large number of 
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13|. The isospin degree of freedom 



enters into the calculations via both cross-sections and mean field 
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14j . The details about 



the elastic and inelastic cross-sections for proton-proton and neutron-neutron collisions can 



be found in Ref. 
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In this model, baryons are represented by Gaussian-shaped density distributions 

1 



Mr,P,t) 



TT 



(1) 



Nucleons are initialized in a sphere with radius R = 1.12 A^^^ fm, in accordance with the 
liquid drop model. Each nucleon occupies a volume of h^, so that phase space is uniformly 
filled. The initial momenta are randomly chosen between and Fermi momentum(pi?). The 
nucleons of target and projectile interact via two and three-body Skyrme forces and Yukawa 
potential. The isospin degree of freedom is treated explicitly by employing a symmetry 
potential and explicit Coulomb forces between protons of colliding target and projectile. 
This helps in achieving correct distribution of protons and neutrons within nucleus. 
The hadrons propagate using Hamilton equations of motion: 

d{ H) dpi d{H ) 
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The baryon-baryon potential V''\ in the above relation, reads as: 

\ J Skyrme ' Yukawa ' C'oul ' ^ sym 

h6{f - f) + t2S{f - f)p^-i (^^^ 
^ ^ exp{\f' - f|//i) ^ ZjZje^ 

3 _ ^^^^ _ 

+t,-nTi6{r^-r^). (4) 

Here and Zj denote the charges of i^^ and j^^ baryon, and Tg, Tg are their respective 
T3 components (i.e. 1/2 for protons and -1/2 for neutrons). Meson potential consists of 

Coulomb interaction only. The parameters fi and ti, ,tQ are adjusted to the real part of 

the nucleonic optical potential. For the density dependence of nucleon optical potential, 
standard Skyrme-type parameterization is employed. The choice of equation of state (or 
compressibility) is still controversial one. Many studies advocate softer matter, whereas, 
much more believe the matter to be harder in nature [12I, [l^ . We shall use both hard (H) 
and soft (S) equations of state that have compressibilities of 380 and 200 MeV, respectively. 

The binary nucleon-nucleon collisions are included by employing the collision term of 
well known VUU-BUU equation {2!, . The binary collisions are done stochastically, in a 
similar way as are done in all transport models. During the propagation, two nucleons are 
supposed to suffer a binary collision if the distance between their centroids 



\ri -rj\ < )J^,<ytot = <y{Vs,type), (5) 

"type" denotes the ingoing collision partners (N-N, N-A, N-tt,..). In addition, Pauli blocking 
(of the final state) of baryons is taken into account by checking the phase space densities 
in the final states. The final phase space fractions Pi and P2 which are already occupied 
by other nucleons are determined for each of the scattering baryons. The collision is then 
blocked with probabihty 

Pblock = 1-(1-Pl)(l-P2). (6) 

The delta decays are checked in an analogous fashion with respect to the phase space of 
the resulting nucleons. 
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III. RESULTS AND DISCUSSION 



The global stopping in heavy-ion collisions has been studied with the help of many dif- 
ferent variables. In earlier studies, one used to relate the rapidity distribution with global 
stopping. The rapidity distribution can be defined as {4, 16 1: 



where E{i) and Pzii) are, respectively, the total energy and longitudinal momentum of i*'* 
particle. For a complete stopping, one expects a single Gaussian shape. Obviously, narrow 
Gaussian indicate better thermalization compared to broader Gaussian. 
The second possibility to probe the degree of stopping is the anisotropy ratio (R) [3]: 

-(EJPii(^)l)' 

where, summation runs over all nucleons. The transverse and longitudinal momenta are 



P±ii) = yPli^) p\\{i) = Pz{i)^ respectively. Naturally, for a complete stopping, 

R should be close to unity. 

Another quantity, which is indicator of nuclear stopping and has been used recently, is the 
quadrupole moment Qzz, defined as [3]: 

Q- = E(2p'W-p'»-pS«)- (9) 

i 

Naturally, for a complete stopping, Qzz should be close to 0. 

In the present analysis, thousands of event were simulated for the neutron-rich reaction 
of 54X6^^^ + 54X6^^^ at incident energies between 50 and 1000 MeV/nucleon using hard 



equation o: 
sections 



state along with energy dependent Cugnon and constant nucleon-nucleon cross- 
171]. Moreover, to see the effect of compressibilities on nuclear stopping and frag- 
mentation, soft equation of state is also used in FigJSl The geometry of the collision was 
varied between the most-central to peripheral one. The role of symmetry energy is studied 
by simulating the above reaction with and without this term. As stated above, we plan to 
study the degree of stopping and emission of fragments using symmetry energy and isospin- 
dependent cross-section. We shall also correlate the degree of stopping with the emission of 
light charged particles as is also done in Ref . ^ . The fragments are constructed within min- 



imum spanning tree (MST) method |l|, which binds nucleons if they are with in a distance 




FIG. 1: The final phase space of a single event for the reaction of 54X6^^^ + 54X6^^^ with(ii) 
and without symmetry energy(i). The top(a), middle(b) and bottom (c) panels are, respectively, 
for scaled impact parameters 6 = 0, 0.3, 0.6. Different symbols are for free nucleons, LCP's and 
IMF's. 

of 4 fm. 

In FigJH we display the final phase space of a single event of 54X6^^^ + 54X6^^^ at inci- 
dent energy of 400 MeV/nucleon, with and without symmetry energy. The top, middle and 
bottom panels are at 6 = 0, 0.3 and 0.6, respectively. Here phase space of free particles [A 
= 1], light charged particles (LCP's) [2<A<4] and intermediate mass fragments(IMF's)[5< 
A < 44] is displayed. We note that irrespective of the symmetry energy, central collisions 
lead to complete spherical distribution of particles, indicating, spreading of the nucleons in 
all directions. It means that breaking of initial correlations among nucleons is maximal in 
this region and, as a result, more randomization and stopping in the hot and compressed 
nuclear matter occurs. This effect seems to decrease with impact parameter. Since free 
particles as well as LCP's originate from the mid-rapidity region, they are better suited for 
studying the degree of stopping reached in a heavy-ion collision. On the other hand, IMF's 
seems to originate either from the target or from the projectile region, therefore, are the 
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FIG. 2: The rapidity distribution ^ as a function of reduced rapidity for free nucleons(a), LCP's 
(b) and IMF's(c) at different impact parameters. The reaction under study is 54X6"'^'^"'^ + 54X6^^^ 
at incident energy E = 400 MeV/nucleon. The left and right panels are with (ii) and without 
symmetry energy (i). 

remnant/residue of the spectator matter. This observation is in agreement with many other 



studies 



18|. 



To further quantify this observation, we display in FigJJl the rapidity distribution ^ for 
the emission of free nucleons as well as LCP's and IMF's. We see that free particles and 
LCP's emitted in the central collisions form a single narrow Gaussian shape, whereas, IMF's 
have broader Gaussian, indicating less thermalization. As we increase impact parameter, 
single Gaussian distribution splits into two Gaussian (at target and projectile rapidities), 
indicating correlated matter. From the shape of the Gaussian, one sees that free particles 
and LCP's are better indicator of thermal source. Obviously, this condition is necessary, 
but, not a sufficient one. 

From the figure, it is also evident that the symmetry energy does not plays significant role 
for the rapidity distribution. The peak value of the Gaussian for LCP's is altered by about 
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FIG. 3: The anisotropy ratio R (i) and quadrupole moment (Qzz/iT'Ucleon) (ii) as a function of 
normalized impact parameter with and without symmetry energy. In panel (b), the results are 
also displayed with isospin-dependent cross-section (55mb). The panels from top to bottom are at 
incident energy of 50 (a), 400 (b), 600 (c) and 1000 MeV/nucleon (d), respectively. 

10%, whereas, nearly no effect is seen in the case of intermediate mass fragments. The 
reason is that LCP's can feel the role of mean field directly, while, the heavy fragments have 
weak sensitivity [Id]. From the figure, one sees a one to one relation between the degree of 
stopping and emission of LCP's. These conclusions match with the findings of FigJT] and 
Ref.|4|. 

In Fig 131 we display impact parameter dependence of global variables {R and Qzz/nucleon), 
whereas, the multiplicity dependence of free nucleons and LCP's is displayed in FigJH The 
displayed results are at Esym= and Esym = 32 MeV in each panel, while, in panel (b) the 
results are also displayed with isospin-dependent cross- sect ion. The value of cross-section 
is denoted in the superscript. From Figj3l We observe that R and Qzz/nudeon behave 
in opposite fashion i.e. R and 7^ — r^— ; — will behave in a similar fashion. For R > 1 and 
Qzz/nudeon < 0, it can be explained by the preponderance of momentum flow perpendicular 
to the beam maximum stopping is observed around 400 MeV/nucleon, 
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FIG. 4: Same as in Figj3l but for the multiplicity of free particles/nucleon and LCP's/nucleon. 



which is in supportive nature with the findings of W. Reisdorf et a/.[2l|. In their work, 
they measured the nuclear stopping from 0.090 to 1.93 GeV/nucleon and maximal stopping 
was observed around 400 MeV/nucleon. It is clear that if the reaction reaches the maximal 
stopping around certain energies, the matter formed in the reaction should reach minimum 
transparency and thus most of the particles are preferentially out-of-plane. On the other 
hand, no visible effect is seen for symmetry energy term. We see both quantities are nearly 
independent of the symmetry energy, while, strongly depends on the isospin-dependent cross- 
section. 

As we know, major contribution for the stopping of nuclear matter is from the hot and 
compressed region, where symmetry energy does not play any role. Some small spikes can 
be seen at lower beam energies, however, outcome is independent of the symmetry energy at 
higher incident energies. This is due to the fact that above the Fermi energy, incident energy 
itself is sufficient to break the initial correlations among the nucleons. On the other hand, 
isospin-dependent cross-section will lead to violent N-N collisions, which further cause the 
transformation of the initial longitudinal motion in other directions and hence thermalization 
of the system. This dominant role played by the isospin-dependent cross-section gradually 
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FIG. 5: Impact parameter dependence of (a) R, (b) Qzz/nudeon, (c) free particles/nucleon, (d) 
LCP's/nucleon, with hard (H) and soft (S) equations of state. The results are displayed at Esym = 
(i) and 32 MeV (ii). 

disappears with increase in the impact parameter. As discussed earher, stopping is the 
phenomena which originates from the participant zone and this zone goes on decreasing with 
increase in the impact parameter and hence the effect of cross-section on nuclear stopping. 
These findings are also in supportive nature with the findings of Liu et al, [3]. 

To correlate the degree of stopping with the multiplicity of fragments, we display in 
FiglU the impact parameter dependence of the multiplicity of free nucleons as well as of 
LCP's. The behavior of all curves is similar to that of nuclear stopping parameters R and 
Q — /J^efeo^ , as discussed in Figj3l In addition, LCP's are more sensitive towards symmetry 
energy compared to free particles. Due to pairing nature of LCP's, symmetry energy term 
oc (A^ ~ Z)"^ contributes considerably. The effect of isosp in-dependent cross-section is more 
visible for the LCP's as compared to free particles. This also gives us clue that LCP's 
production can act as a indicator for the nuclear stopping. Moreover, free particles/nucleon 
are found to increase monotonically with the incident energy, while LCP's/nucleon behave 
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FIG. 6: System size dependence of (a) R, (b) g — /nudeon ' ("^) ^^^^ particles, (d) LCP's, in the 
presence of symmetry energy. All the curves are fitted with power law. 

in similar fashion as that of nuclear stopping i.e. maximum around 400 MeV/nucIeon and 



then decreases. It is also evident from Ref.j4|, LCP's production act as a barometer for 
nuclear stopping compared to the free particles. 

In Figini we have checked the sensitivity of nuclear stopping as well as fragment pro- 
duction with the nuclear equation of state (EOS). For this purpose, a hard (H) and soft 
(S) equations of state with compressibility k = 380, 200 MeV are employed, respectively. 
The nuclear stopping is found to be weakly dependent on the equations of state, while, the 
fragments production is sensitive to different equations of state. It means that the fragment 
production with different equations of state can act as a global indicator for the nuclear 
stopping as it is weakly dependent on equations of state. 

It also becomes important to study the system size dependence and isospin asymmetry of 
R, Q /nudeon ' ^^^^ particlcs and LCP's. For this, in FigJHl we have displayed the results for 
the reactions of 2oCa^° + 2oCa^°, 2&Nr'^ + 2&Nr'^, ^iNb^^ + 4lA^&^^ 54Xe^^^ + ^iXe^^^ 
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FIG. 7: Isospin asymmetry of (a) R, (b) -q — /nudeon ' ^^^^ particles, (d) LCP's, in the presence 
of symmetry energy and isospin-dependent cross-section (55 mb). 

and 79^4^^^''' + jqAu^^"^ , in which Z as well as A is varied. On the other hand, results are 
displayed, in FigJH for the reactions of aoCa^^ + aoCa^^ (N/Z = 0.7), 2oCa^° + 2oCa^° 
(N/Z = 1), 2oCa^^ + 2oCa^® (N/Z = 1.4) and 2oCa" + 2oCa^^ (N/Z = 1.85), having same 
Z and different A, in the presence of symmetry energy and isospin-dependent cross-section. 
The curves in Figsl6] and [7] are parametrized with the power law Y = CX'^, where C and r 
are constants, while X and Y are the respective parameters on X and Y axis. 

From the Fig|6l it is observed that the parameters R, g — /nucieon ' ^^^^ particles as well 
as LCP's are in similar trend with the composite mass of the system. All the parameters 
are found to increase with the composite mass of the system. For a fixed geometry (semi- 
central here), more heavier is the composite system, more hot is the compressed zone, which 
further results in more thermalization or global stopping. Looking the parallel side, the free 
particles and LCP's will always originate from the participant zone. With an increase in the 
composite mass of the system, the participant zone goes on increasing for a fixed geometry 
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FIG. 8: The scaled free particles/nucleon and LCP's/nucleon along with anisotropy ratio R (a) and 
Q — jnucleon ^ ^ function of normalized impact parameter in the presence of symmetry energy. 
The reaction is at incident energy 400 MeV/nucleon. 



(semi-central here) and hence the production of free particles and LCP,s. Similar findings 
are also published in the Ref.[3, 122I . 

The dependence of these parameters on the isospin asymmetry (N/Z dependence ) displayed, 
in FigJTl is also found to be in supportive nature with the findings in FigJHl An increase in 
the number of neutrons will increase the number of collisions and hence dominance of R, 
7^ — r^-^ — , free particles as well as LCP's is observed with increase in N/Z ratio. Nuclear 
stopping as well as LCP's are observed to be strongly dependent on the isospin-dependent 
cross-section. Similar results with isospin dependent cross-section are observed in Figsj3]and 
m From here, one may conclude that the nuclear stopping and LCP's can also be used as a 
tool to investigate the isospin-dependent cross- sect ion. 

To further elaborate this point, we display in FigJHl Multiplicity/nucleon (free and 
LCP's) as well as R and 7=; — r^— ; — . Once free nucleons and LCP's are normalized with 

' Qzz/nuclecm 
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R at the starting point of impact parameter, we see that their behavior with respect to 
impact parameter is similar to that of anisotropy ratio, whereas, visible difference occurs 
with reference to quadrupole moment. This similarity in all three quantities in the presence 
of symmetry energy makes LCP's good indicator of global stopping in heavy-ion collisions. 

IV. CONCLUSION 

In summary, using the isospin-dependent quantum molecular dynamics (IQMD)model, 
we investigate the emission of free particles, LCP's, and degree of stopping reached in a 
heavy-ion collisions in the presence of symmetry energy and isospin dependent cross-section. 
We observed that nuclear stopping in term of anisotropy ratio and quadrupole moment de- 
pends weakly on the symmetry energy and strongly on the isospin-dependent cross-section. 
On the other hand, the symmetry energy and isospin-dependent cross-section has an effect 
of 10% on the LCP's production. It means nuclear stopping and LCP's production can be 
used as a tool to investigate the isospin-dependent cross-section. The LCP's production is 
found to be highly correlated with the global stopping. 
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